Effective spin-flip scattering in diffusive superconducting proximity systems 

with magnetic disorder 



D. A. Ivanov,^ Ya. V. Fominov,-^ M. A. Skvortsov,^ and P. M. Ostrovsky^'^ 

^Institute of Theoretical Physics, Ecole Polytechnique Federale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 
D. Landau Institute for Theoretical Physics RAS, 119334 Moscow, Russia 
'^Institut fiir Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 

(Dated: 5 October 2009) 

We revisit the problem of diffusive proximity systems involving superconductors and normal met- 
als (or ferromagnets) with magnetic disorder. On the length scales much larger than its correlation 
length, the effect of sufficiently weak magnetic disorder may be incorporated as a local spin-ffip term 
in the Usadel equations. We derive this spin-flip term in the general case of a three-dimensional 
disordered Zeeman-type field with an arbitrary correlation length. Three different regimes may be 
distinguished; pointlike impurities (the correlation length is shorter than the Fermi wavelength), 
medium-range disorder (the correlation length between the Fermi wavelength and the mean free 
path), and long-range disorder (the correlation length longer than the mean free path). We discuss 
the relations between these three regimes by using the three overlapping approaches: the Usadel 
equations, the non-linear sigma model, and the diagrammatic expansion. The expressions for the 
spin-fiip rate agree with the existing results obtained in less general situations. 

Published as: Phys. Rev. B 80, 134501 (2009) 

PACS numbers: 74.45.-l-c, 75.60.Ch, 74.78.Fk 



Contents 



I. Introduction U 

II. Main results [2 

A. Parameters of the problem 

B. Spin- flip term in the sigma model 

C. Spin- flip term in the Usadel equations 

D. Spin-flip scattering rate 

1. Regime of short-range correlations 

2. Regime of medium-range correlations 

3. Regime of long-range correlations 

III. Sigma-model derivation 

A. Sigma-model action 

B. Local spin-flip term 

C. Nonlocal spin- flip term 

D. Discussion of the two contributions 

IV. Derivation from the Usadel equations 

A. Effective spin-flip rate 

B. Self-consistent screening of disorder- induced 
correlations Q 

C. Effective spin- flip scattering at the edge of the 
minigap 

V. Diagrammatic representation 

VI. Conclusions [l^ 

Acknowledgments [l^ 

References [l^ 



I. INTRODUCTION 

In conventional superconductors (S) , pairing occurs be- 
tween electrons with opposite spins, and thus, the co- 
existence of superconductivity and magnetism may lead 
to a variety of interesting effects in superconducting 
and proximity structures. Examples of such effects are 
gapless superconductivity,-i triplet proximity correlations 
(see Ref. for a review), and Josephson 7r-junctions (see 
Refs. for a review). In a ferromagnet (F) with a uni- 
form exchange field, the theory of anomalous correlations 
may be constructed by taking into account the splitting 
of electronic bands of opposite spins*2. The situation be- 
comes more complicated if the magnetic structure is in- 
homogeneous. In this case, if the ferromagnetic structure 
is non-coUinear, the triplet component of anomalous cor- 
relations needs to be taken into consideration!^ The case 
of inhomogeneous magnetic structure is also interesting 
from the practical point of view, since many experimen- 
tal studies of hybrid superconductor-ferromagnet (SF) 
systems, in particular, Josephson 7r-junctions, reveal a 
strong spin-flip scatteringi^iii^ A theoretical analysis of 
superconducting correlations in the presence of an inho- 
mogeneous magnetic structure is complicated, and many 
of the existing studies are limited to considering either 

pointlike magnetic impurities^'^'^ i^^'^^i^^'^'^ or specific do- 
main geometries^^^^^^^^^^^^^^26,27 

In a recent work by two of the present authors, 
the spin-flip term (in the same form as for magnetic 
impurities^) has been derived in the model of non- 
collinear magnetic disorder correlated on length scales 
much larger than the elastic mean free pathJS The re- 
sulting expression for the spin-flip rate agrees with the 
known one obtained for the coUinear inhomogeneous 
magnetization under assumption of periodic magnetic 
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FIG. 1: Two types of spin-flip scattering: (a) local; (b) non- 
local. The diagrams show the way the magnetic scattering is 
included in the "cooperon" propagator (see Sees. IIIII and |V] 
for details). Thick solid lines represent the electron Green 
functions, thin dashed lines denote averaging over Gaussian 
disorder. The crosses and the open circles represent the po- 
tential and magnetic disorders, respectively. 



well as the expressions for the spin-flip scattering rate 
in the three regimes, including the crossovers between 
those regimes. In the following sections, we present de- 
tails of the calculations. In Sec. IIIII we employ the sigma 
model to study all cases of the magnetic-disorder correla- 
tion length except for the crossover between the medium- 
and long-range regimes. In Sec. IIVI we demonstrate how 
the long-range regime can be treated in the language of 
the Usadel equations. In Sec. |Vl we consider all the cases 
with the help of the diagrammatic technique. Finally, in 
Sec. IVII we present our conclusions. 



struct ureiiiiii 

In the present paper, we revisit the problem of the 
magnetic disorder in the more general situation of the 
non-coUinear disorder with arbitrary correlation length. 
We consider a superconducting or proximity-type system 
with potential impurities and an inhomogeneous Zeeman 
field. The potential impurities are supposed to be suffi- 
ciently strong to bring the electronic motion to the dif- 
fusive regime. On top of this diffusive motion, the elec- 
trons experience splitting from the inhomogeneous Zee- 
man field, which is assumed to be random and Gaussian 
with an arbitrary pair correlation. We further assume 
that this magnetic disorder is much weaker than the po- 
tential one, in terms of the characteristic scattering rates. 

Then three different regimes can be distinguished: the 
short-range magnetic disorder (or, equivalently, point- 
like impurities, with the correlation length shorter than 
the Fermi wavelength), the medium-range disorder (with 
the correlation length between the Fermi wavelength and 
the elastic mean free path), and the long-range disorder 
(the correlation length longer than the elastic mean free 
path). The short-range case has been solved in Ref. ilj, 
the medium- and long-range regimes have been treated 
m Refs. [Hill for the coUinear periodic case, and the 
long-range non-coUinear case was studied in Ref. [2^. We 
extend those results to the general non-coUinear case 
and remove some of the technical assumptions made in 
Ref. M 

We use three methods for our analysis: the non-linear 
sigma model, the Usadel equations, and the diagram- 
matic technique. While the calculations in these three 
methods are somewhat parallel to each other, we find 
it instructive to present those various approaches, in or- 
der to illustrate the correspondence between the meth- 
ods and to clarify the physical meaning and the appli- 
cability conditions of the results. As we shall see below, 
the spin-flip term in the short-range and medium-range 
regimes corresponds to inserting one magnetic-impurity 
rung into the cooperon ladder [Fig.[Ha)], while the long- 
range regime corresponds to the magnetic line crossing 
many cooperon rungs [Fig.[TJb)]. 

The paper is organized as follows. In Sec.|TTl we present 
the main results of the paper: the form of the spin- 
flip term in the sigma-model and Usadel description, as 



II. MAIN RESULTS 
A. Parameters of the problem 

We consider a diffusive motion of electrons in a finite 
piece of magnetic metal (i.e., a metal with exchange in- 
teraction inside of it) of linear size L. The Thouless en- 
ergy scale (the inverse diffusion time) of this region is 
Etu = D/L'^, where D is the diffusion constant (we put 
h = 1). Below we study the anomalous correlations in 
this magnetic metal induced either by a small (see condi- 
tions below) superconducting order parameter A or by 
an electric contact with a superconductor (proximity- 
induced superconductivity, in which case the order pa- 
rameter A is put to zero in the magnetic metal). The 
electrons and the (Andreev-reflected) holes are consid- 
ered at a finite energy E (relative to the Fermi level). 
We further consider an exchange field in the magnetic 
metal H(r) = h(r) + Sh{r) containing a small (possibly 
zero) smooth (varying on the length scale L) part h and 
the disorder component Sh (the typical scale of the dis- 
order part will be further denoted as 5h). The magnetic 
disorder is assumed to be correlated on a short length 
scale a, which defines the scale of "the Thouless energy 
of magnetic inhoniogeneities" Ea — D/a?. We assume 
a Gaussian ensemble for (5h, with the pair correlation 
function 

{5h,{v)5h,{v'))^F,,{\v-v'\). (1) 

The typical order of magnitude for Fij{r) is then (Sh)'^, 
and the typical support is of order a. 

For the energy scales in the magnetic metal, the fol- 
lowing condition is assumed)^ 

ETh,E,A,h,T,f^Ea,T-\ (2) 

where Fgf is the resulting effective spin-flip rate, and t is 
the mean free time due to potential scattering. The phys- 
ical meaning of this condition is that the length scales as- 
sociated with both potential and magnetic disorder (the 
right-hand side of the inequality) are much shorter than 
the length scales involved in the Usadel equations (the 
left-hand side of the inequality). 
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B. Spin-flip term in the sigma model 

The sigma-model action has the form 

S[Q]=So+Ssi, (3) 
where the usual sigma-model action is 

+ (zSfg- A-ih(r)f36-)g} (4) 

(we use the standard notations of the sigma-model tech- 
nique, see Sec. Illll for definitions), and the spin-flip term 
is 



j d?v T% STr (fga.Qfaa.Q) , (5) 



where F^j is a symmetric matrix of the spin-flip scatter- 
ing rates Here and below, we use the convention of 
summing over repeating indices. 



terms: 
D 



V'^e + Mo (iE sin 61 -I- A cos 6) - (hM) cos 9 
- (rlf + 2M'^f m) sin 29 = 0, (9) 



y (MV^A-fo - A-foV^M) - M(i£:cos6' - Asin6i) 

- Moh sin 9 + 2A/of sf M cos 29 = 0, (10) 

where Tsf is the symmetric 3x3 matrix of F*], and 
ptot _ TrFsf. This set of equations generalizes the con- 
ventional spin-flip termSi to the case of the triplet Usadel 



equations 



30 



D. Spin-flip scattering rate 

1. Regime of short-range correlations 
In the regime of a <C A:^^, the spin-flip rates are given 



F^^f = TTi/y d^r Fij{r) - {5hf . 



(11) 



C. Spin-flip term in the Usadel equations 

The spin-flip term in the Usadel equations can be ob- 
tained by varying the action ^ with respect to the Q 
matrix. Denoting the saddle-point value of the Q matrix 
as the matrix Green function g, we write the resulting 
equation as 



-DV {gVg) + iEf^ao - Acto - if3(h<T), g 



[Tso-jgfao-j ,^] = 0. (6) 



The non-linear constraint g^ = \ can be resolved by 
a parametrization. In the case when the phase of the 
anomalous Green function is fixed (in which case the Us- 
adel equations contain four parameters), we may use the 
parametrization of Ref. i28| in terms of the angle 9 and 
the vector M: 



g — Mf)(Ta (fs cos 9 + fi sin 9) + iMo-(f3 sin 9 — fi cos ( 



with the constraint 



- = 1. 



(7) 



(8) 



A generalization to the case of a varying phase is straight- 
forward, in which case four more parameters need to be 
added, corresponding to the terms containing f2 in the g 
matrix. 

Substituting the parametrization ([7]) into Eq. ([5]), we 
obtain the triplet Usadel equations with the spin-fiip 



2. Regime of medium-range correlations 

In the regime oi kp^ a ^ I, the spin-flip rates are 
given by 
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F^^TTi. J d'rF,,{r)^^^^^{Shfakf. (12) 

The exact formula interpolating between the two 
regimes ([TT]) and has the form 



rl^, = n,y d'vF,,{r) 



'3„ I ^sin^(fci=^0 



(13) 



3. Regime of long-range correlations 
In the regime of Z <C a, the spin-flip rates are given by 

-,2 



r- = ^JdhF,,ir)^^^.iShf 



I 



(14) 



The exact formula interpolating between the two 
regimes (IT^ and has the form 



dh 



3D J (27r)3 



arctan(g/) 
ql — arctan(g/) 



(15) 



Note that the diffusion constant is related to the den- 
sity of states by £> = vpl/3 = kpl/{6TT^i'). 

The results ([T4|) and have been derived for 

collinear periodic magnetic structures in Refs. fl^ l Tsl The 
result P^ in the non-collinear isotropic case has also 
been found in Ref. u^. 
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III. SIGMA-MODEL DERIVATION 

In the sigma-model description, if we assume that the 
random exchange field Sh is sufficiently weak, then it 
can be included in the sigma-model action perturbatively. 
For the delta-correlated (short-range) magnetic disorder, 
this procedure is well known ; and our consideration 
generalizes it to arbitrary correlation lengths. 

In our derivation, we find that there are two different 
contributions arising from the magnetic disorder. By ex- 
panding the action in the disordered field to the second 
order. 



S — Sq + Si + S2 ■, 
we find the two spin-flip-type contributions: 



(16) 



1 



^loc = {S2) and 5„o„ioc = -2^^'^^''^) ' ^'^'^^ 

We shall see below, that the first ("local") contribution 
dominates in the case of short-range or medium-range 
disorder and gives the spin-flip rate pS]) . while the sec- 
ond ("nonlocal") contribution becomes dominant in the 
regime of long-range disorder and gives the spin-flip rate 
(I14p . In performing the averaging for the nonlocal con- 
tribution, one needs to take into account fluctuations 
around the replica-symmetric (or supersymmetric) sad- 
dle point. On inspection, the local and nonlocal spin- 
flip contributions correspond to the processes depicted 
in Fig. [1] (left and right panel, respectively). Interpolat- 
ing between these two regimes goes beyond the scope of 
the sigma-model derivation in this Section, but it is done 
in Sec. |V]in the diagrammatic language. 



A. Sigma-model action 

To derive the sigma model for a disordered super- 
conducting (or proximity) system, we follow the usual 
procedureJ^i^ The potential disorder is assumed to be 
Gaussian and 5-correlated, 



(C/(r)[/(r'))j 



^(r - r') 

2'KVT 



(18) 



[here v = rnkp / {'2.tt'^) is the density of states at the Fermi 
level per one spin projection], while the magnetic disor- 
der is also taken to be Gaussian, but with an arbitrary 
correlation length, see Eq. ([1]). 

As usual in the derivation of the sigma model, we as- 
sume the "dirty limit", i.e., that 



£;Th,-B,A,;i<T- 



(19) 



In addition, we need that the effective spin-flip scattering 
(whose rate we derive below) is much weaker than the 
potential scattering. 



Fsf < T- 



(20) 



so that we can treat it as a perturbation on top of the 
diffusive sigma model defined by the potential disorder. 

The derivation of the sigma-model action starts with 
the partition function for excitations with energy E in 
the Bogolyubov - de Gennes Hamiltonian including also 
the exchange field)^ 



Z = / D^*D^ e 



-S 



<S[*] = -i j 



E~na-TJi + U{v) 



f2 Re A - fi Im A 



2m dr^ 



(21) 



(22) 



(23) 



Here, ^* and ^' are 4-component (in the product of 
the Nambu-Gor'kov and spin spaces) fermionic vector 
fields containing Grassmann anticommuting elements. 
For brevity, we do not write a small imaginary part iO 
that should be added to the energy E. The Pauli ma- 
trices ai act on the spin of electrons while fi act in the 
Nambu-Gor'kov space. 

Now we introduce replicas^^^ (or supersymmetryi^) 
for averaging over the potential disorder, thus extending 
the vector. As a result of averaging with the disorder 
correlator ITSl). we obtain 



(fir 



fi Im A 



£; - Ho- - faC - f2 Re A 
* + / d'r (*+f3*) 



(24) 



Next, we perform the Hubbard-Stratonovich trans- 
formation with the help of a matrix field Q{r). The 
Q field is a matrix in the product of the Nambu- 
Gor'kov and spin spaces and also in the space of repli- 
cas (or in the Fermi-Bose superspace). As a result of 
the Hubbard-Stratonovich transformation, the action be- 
comes quadratic, and the Gaussian integration over 4"* 
and ^ yields 



S[Q] = ^ j <fv STrg2 



<Fv STrln 



e - f3 (i? - H^) - *A - g 



A = 



A 

A* 



, (25) 
(26) 



The "STr" operation is the complete trace. In particu- 
lar, it includes the trace over the replica indices or the 
supertrace, depending on the version of the sigma model 
(replica or supersymmetric). 

In general, the Q matrix of a sigma model for a su- 
perconducting system would contain all the diffusive soft 
modes: diffusons and cooperons. This would be achieved 
by an additional doubling of the fermionic fields /^'^'^i'^^ 
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and hence of the Q matrix, in the retarded-advanced 
space. The doubled Q matrix would be subject to an 
additional constraint reflecting the symmetry of the Bo- 
golyubov - de Gennes Hamiltonian. Note that our sigma 
model does not involve this doubling and accordingly in- 
cludes only the cooperon modes but not the diffusons. 
As a result, our calculations are valid only at the saddle- 
point level (Usadel equations), but cannot reproduce 
weak-localization corrections involving diffuson modes. 
An interesting question of the effect of the inhomoge- 
neous magnetic disorder on the weak-localization correc- 
tions requires a more detailed analysis and is postponed 
until further studies. In this paper, since we are inter- 
ested in the effect of the magnetic disorder on the saddle 
point, we use the reduced version of the sigma model. 

In the quasiclassical regime, where the Fermi energy 
is the largest energy scale, Epr 3> 1, the Q matrix is 
restricted to the manifol d^^'^^ 



B. Local spin- flip term 

Averaging S2 over magnetic disorder ([1]) produces the 
spin-flip term ^ with 

T^^..j^rF.,ir^^e-^n, (33) 

Note that this expression generalizes the result of 
Abrikosov and Gor'kov for pointlike impurities.™ In 
that work, the magnetic disorder was assumed delta- 
correlated, which corresponds to ^^(r — r') = 5ij5{Y — 
y') / {QiTVTs) and = (5y /(6Ts). Our derivation extends 
that result to the medium-ranged disorder with correla- 
tion lengths up to L As we shall see below, the contribu- 
tion (j33p is dominant as long as a ^ /, and therefore in 
this regime we can neglect the factor e^*^^' and arrive at 
Eq. USD. 



(27) 



Furthermore, using the dirty- limit assumption (|19p . we 
expand the action in the gradients of Q [simultaneously 
expanding the logarithm in Eq. (pS)) in 5\\\ and obtain, 
in the usual manner i^^i'^^ the action (flHt with 



Sq=ttv / d^r STr { ^{VQ f + ih{E - ^h) - A Q 



Si — —inu 



J d^rdhir) STr(f36-g(r)), 



(28) 
(29) 



and 



5, = - STr / d-^rrf^r' 



fd^pd^p' 
(27r)6 



[TT3&5hT-^] 



gip(r-r') 
C iA 

g-ip'(r-r') 

x[Tf,aShT-']^, J, , (30) 

? 2r 



where the local matrix T(r) parametrizes rotations of the 
Q matrix, 

g^T-^AT, A = f3, (31) 

and ^ = pV(2m) - = P'V(2"^) - 

The integrals over p and p' in Eq. ([50]) may be com- 
puted, giving rise to the kernel decaying at the elas- 
tic scattering length I. Assuming that the Q matrix 
changes on length scales much longer than ?, we can put 
T(r) = r(r') in Eq. ^ and arrive at 



S2 



2 J {kF\r-r'\Y 
X 6h^{r)6hj{r') STr(f3<7.Q(r)f3(7jQ(r)) . (32) 



C. Nonlocal spin-flip term 

Averaging Si over the magnetic disorder ([T]) produces 
the contribution to the action 

^^((5i)^) = ^/dWr%(r-r') 

X STr(f3<T,Q(r)) STr(f3ajQ(r')) . (34) 

While the main part of the action Sq contains only 
one STr, this contribution is a product of two super- 
traces. We assume that the saddle point Qq is super- 
symmetric (or replica symmetric), and then the con- 
tribution of Eq. vanishes at such a saddle point. 
However, taking into account non-supersymmetric (non- 
replica-symmetric) fluctuations around the saddle point 
produces a non-negligible contribution containing only 
one STr. 

In order to average Eq. ([34]) over fluctuations of Q, we 
parametrize those fluctuations by local rotation matrices 
W , anticommuting with Qq: 



Q^Qo + iQoW 



(35) 



The effective action for W, extracted from the Sq part, 
to the Gaussian order is 



■nvD 



j d^rSTviVWf 



(36) 



Note that since Eq. (l34|) involves correlations of W at 
the length scale of order a, we only need to take into 
account short-wavelength fluctuations of W . Therefore, 
we neglect the terms containing h. A, and VQo in 
Eq. ([36|) . as well as the self-consistent "screening" by the 
effective spin flip Est [which produces the infrared cutoff 
for the action ([55]) : see also Sec. lIVBj below]. under the 
assumption 



(37) 



To average over the fluctuations with the action pH)) 
for W anticommuting with Qq, we employ the following 
contraction rule>i^'^ 

(STr(AiW^(r)) STr(A2Vt^(r')))5^^, 

= -^[V-\r' STt{AiQoA2Qo - A1A2) (38) 

for any operators Ai and A2 [here [V^^Jrr' ~ — (47r|r — 
r'|)~^ is the kernel of the inverse Laplacian in three di- 
mensions]. Applying this identity to averaging the con- 
tribution (|34p with the matrix Q parametrized by Eq. 
([35)) . one finds that — i ((iSi)^)^ is given by the usual 
spin-flip term ([5]) with 

^7f = -^ I rf'rF,, (r)[V-V, (39) 
which coincides with Eq. (fT4|) . 

D. Discussion of the two contributions 

While the two contributions fTT]) come from two dif- 
ferent terms in the sigma-model action, we shall see that 
they, in fact, correspond to the two limiting cases of mag- 
netic disorder, as depicted in Fig. [1] 

First of all, let us compare the magnitude of the two 
spin-flip rates. The "local" spin-flip rate ([33|) has the 
order of magnitude 

r'^°''''^ ^i^iShfa', a^kp\ (40) 
r^;°^'^') ^ 1^ {5hf akf , fc^i « a < / , (41) 
T^^r^'^^'^i^hflkf, Z«a. (42) 

On the other hand, the "nonlocal" contribution ([39]) can 
be estimated as 

^(nonlocal) _ _ ^^Shf^kf . (43) 

Therefore, the nonlocal term ([39]) dominates for a ^ I, 
while the local term (|33p becomes dominant at a <ti I 
[strictly speaking, the nonlocal term is only defined for 
I, see our discussion below]. 
Second, we can graphically represent the two contribu- 
tions as shown in Fig. [2] [panels (a) and (b)] . Correlation 
functions of the Q matrices in the sigma model corre- 
spond to the diffusion ladders in the conventional dia- 
grammatic techniquCfi^ and therefore the sigma-model 
diagrams represented in Figs. [2][a) and [U^b) translate 
to the processes shown in Fig. [TJ Thus, while for- 
mally our derivation produces the sum of the two terms 
p(iocai) _|_ p(noniocai) ^ ^^^^ ^j^^ local term should be kept 

in the regime a <^ I, and only the nonlocal term — in 
the opposite regime a ^ I. At the same time, the nonlo- 
cal term has been derived under an implicit assumption 
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(a) (b) (c) 

FIG. 2: Diagrammatic representations of the (a) local and 
(b) nonlocal contributions in the sigma-model calculation. 
Thick double lines correspond to the Q matrix; closed loops of 
thin solid lines indicate the STr operation; open circles denote 
the "magnetic vertices" faff. Dashed lines are the magnetic 
disorder correlation functions Fij; wavy lines are the propa- 
gators of the W field. Panel (c) shows a potentially dangerous 
contribution to the propagator of the W field. To avoid di- 
vergencies, a self-consistent "screening" of the W propagator 
needs to be taken into account. 

a ^ I, since in Eq. psp we have used the diffusion prop- 
agator for the correlation function of W{r) and W{r') 
at distances of order a. Therefore, at the intermediate 
length scales a ^ I, none of the terms ([55)1 or ([55)l . nor 
their sum, provide an accurate result for the spin-flip 
rate. To calculate an effective spin-flip rate at a ~ a 
crossover from the ballistic to the diffusive motion needs 
to be taken into account. This calculation is performed 
in Sec.|V|in the diagrammatic language (a similar calcu- 
lation in the context of a collinear periodic magnetization 
has been done in Refs. fl^lisl ). 

Finally, we would like to comment on the applicability 
conditions of our derivation. For the derivation of the 
local term (at a ^ Z), we only need the conditions ([TO]) 
and ([^0)1 . For the nonlocal term (at a ^ /), we have also 
assumed the condition (|37p . Altogether, the applicability 
conditions may be reformulated in the universal form ([2]). 

Note that the propagator of the W field (used in the 
calculation of the nonlocal term) gets, in principle, renor- 
malized by higher-order contributions. One of the po- 
tentially dangerous corrections is shown in Fig. ^c) . To 
avoid infrared divergencies in this correction, one needs 
to take into account that the diffusive propagator p6]) 
for the W field is cut off at large distances by a certain 
screening length. This screening is discussed in more de- 
tail (in the language of the Usadel equations) in Sec. lIVBl 
the resulting screening length is given by Eq. ([5S]) . If this 
screening is taken into account, then higher-order correc- 
tions to the propagator of the W field may be neglected. 

IV. DERIVATION FROM THE USADEL 
EQUATIONS 

In this section, we present an alternative derivation 
of the spin-flip term in the regime of long-range cor- 
relations (a 3> by directly averaging the Usadel 
equationsSii^^i^i^ over the magnetic disorder, following 
an approach similar to that of Ref. [2^. In this way, we 
derive Eqs. ([6]) and ([T4]) and lift the assumptions of "self- 
averaging" and of "being away from the gap edge" im- 
posed in Ref. |28|. 
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In the regime of long-range correlations, the general 
assumption ^ may be simplified as 

ETi„E,A,h,5h<^Ea (44) 

[since Fsf ^ {Sh)'^/Ea, as we derive below]. 



A. Effective spin-flip rate 

We start from the Usadel equation containing the ex- 
change fieldf^i^ 

DV (GVG) + \iET3a0 - A(7o - iT3(Ua),G] = 0, (45) 



where the gradient term can also be rewritten as 
V(GVG) - i [G,V^G] 



(46) 



due to the normalization condition & = 1. Here, H is 
the total realization-dependent exchange field, containing 
a smooth background field h and a Gaussian disorder Sh 
obeying Eq. ([T]): 

H = h + (5h. (47) 
The exact solution of Eq. (|45p can be written as the sum 



G = g + Sg 



(48) 



of the disorder-averaged part g = (G) and the Sg part 
that depends on the realization of the magnetic disorder 
and averages to zero. 

As confirmed by our further derivation, under the as- 
sumption (|44p , the realization-dependent part Sg is small, 
\Sg\ <^ I, and it is sufhcicnt to consider it linear in Sh. 
Our aim is to obtain the equation for g, the disorder- 
averaged Green function. Note that g is not simply the 
zeroth order over 6h: under our assumptions it is also 
influenced by the averages containing the second order 
over i5h. 

Averaging the normalization condition G^ = 1 over 
the magnetic disorder and neglecting the {Sg"^) term [this 
is possible under the assumption (|44)) . see Sec. IIVB] for 
details], we obtain the normalization condition g^ = 1. 
Then the realization-dependent part must obey the rela- 
tion 



{9,Sg} = 0. 



(49) 



Averaging Eq. (I45p over the magnetic disorder and tak- 
ing into account Eq. (|46|) . we find 



D 



iETs&o - Ao-Q - iT3{h.a-),g 

-i[T3a„{6h,6g)]^0, (50) 



where the summation over the repeating indices is as- 
sumed. Here, we have dropped out the full derivative 



term containing V {[5g, ^5g]), since this term has an ad- 
ditional smallness [as confirmed by the result ([5^ below] . 

To calculate the averages in Eq. ([50| , we extract from 
Eq. (|^5|) [we also take into account Eq. ([^5]) ] the linear 
part in 6h.: 



D 



D. 



= t[MSh&),g]. (51) 



This is a linear equation with respect to Sg, with the 
source term containing the disorder Sh. In order to find 
Sg from this equation, we note that the first term is the 
largest one, since the derivatives apply to the fast func- 
tion Sg which follows Sh and hence changes on the scale 
of a. The second term in the left-hand side is smaller, 
according to our assumption ()44|1 . and for most purposes 
we may neglect it. Employing Eq. ()49|) . we then obtain 



(52) 



^ff = (V '^Shi) g[T3a,,g] 



Now the disorder-induced part of Eq. ([50]) after simple 
algebraic manipulations takes the standard form of the 
spin-flip term: 

- i [Ts&i, (ShiSg)] = -Til [^3^^3x30-^,5] , (53) 

where we have defined 



= --{sh,y-'sh,) 



(54) 



[which is equivalent to Eq. (|14p ]. Thus, we finally arrive 
at Eq. dll). 



B. Self-consistent screening of disorder-induced 
correlations 



In the above calculation, we have neglected the terms 
containing the disorder averages {Sg^), at the same time 
keeping the terms with {ShiSg). While the exact calcu- 
lation of the neglected terms appears to be a delicate 
problem, we can estimate their order of magnitude (and 
thus justify neglecting them) from simple arguments. 

To estimate those averages, we should express Sg from 
Eq. ([5T|) and then average over Sh. However, in this pro- 
cedure, we cannot limit ourselves to the approximation 
(|52p . which would produce a divergence: 



(55) 



To regularize this infrared divergence, one needs to 
take into account the second term in Eq. (|5ip , which pro- 
vides an effective cut-off for the integral ((55|) . The cor- 
responding "screening length" R* is determined by the 
largest of the energy scales E'Th, E, A, h [corresponding 
to the second term in Eq. ([?T|) ]. and Egf [the latter energy 



8 



scale appears if one self-consistently includes the spin-flip 
terms in Eq. ^ in our expansion]: 

^* ^ \/max(i?Th,i?,AXr:0 ■ ^^^^ 
As a result, the neglected terms may be estimated as 

For our approximation (neglecting those terms), we 
require that they are much smaller than one, and that 
they produce corrections smaller than the spin-flip term 
Tgf ^ {5hY /Ea- The first condition translates into the 
requirement 5h <C Ea (under this condition, the ne- 
glected terms are small). The second condition gives the 
additional constraint Et^jE, A,h Ea (under this con- 
dition, the spin-flip term is the main effect of the disor- 
der). Altogether, the conditions of applicability of our 
derivation can now be formulated as Eq. (j44p . 

Note that taking into account the cut-off length R* al- 
lows us to get rid of the assumption of the "self-averaging 
disorder" made in Ref. 5^ to guarantee the convergence 
of the integral (|55| . Now we see that the condition ([44|l 
is sufficient for that. 

The above treatment of the self-consistent screening of 
disorder-induced correlations has been performed under 
the assumption of a three-dimensional magnetic disorder. 
At lower dimensions (e.g., if the magnetic disorder forms 
layers), a nonlocal interference becomes important, and 
replacing the effect of magnetic disorder by a local spin 
flip may not always be possible.— Technically, the three- 
dimensionality is used in our derivation of the estimate 
([57)1 . If one repeats the calculation in a dimension lower 
than three, one finds that the smallness of {Sg'^) cannot 
be guaranteed. So in low dimensions our approximation 
breaks down, and we expect that the role of magnetic 
disorder becomes in this case nonlocal and non-universal. 



C. Effective spin- flip scattering at the edge of the 
minigap 

In the above derivation, we assumed that the linear 
operator acting on Sg in Eq. (j5ip is invertible [and that 
we can keep only the term in its inverse, which led 
us to Eq. ([52]) ] . As it was pointed out in Ref. [23, this as- 
sumption breaks down at the edge of the minigap, where 
the solution of the Usadel equation bifurcates. In that 
case, the linear-order perturbation theory over i5h breaks 
down, which is formally reflected in the non-invertibility 
of the linear operator in Eq. ([5T[) . 

However, the non-invertibility of the operator (and the 
corresponding zero mode) is associated with the bound- 
ary conditions at the edge of the system and, hence, with 
the length scale L of the system size. On the other hand, 
the spin-flip processes producing the scattering rate (fH[) 



and corresponding to Fig. [TJb) are associated with the 
length scale a (with a ^ L) and do not depend on the 
boundary conditions. Therefore we expect that the same 
form of the spin-flip term remains valid also near the 
minigap edge. 

The problem with the perturbative expansion at the 
minigap edge is due to the fact that the minigap itself 
depends on the spin-flip rate,— which leads to a non- 
analytic shift of the saddle point of the action ([3]) as a 
function of the disorder 5\\. A correct way to include 
the effect of disorder in the vicinity of the minigap is to 
first calculate the spin-flip rate from Eq. ([T4[) [this ex- 
pression is independent of the particular saddle-point so- 
lution ^(r)] and then to recalculate the new saddle-point 
solution (with a new value of the minigap) with this spin- 
flip rate using the Usadel equations. We thus conclude 
that our result for the effective spin flip remains valid at 
all ener^s across the minigap edge, and the assumption 
of Ref. 28 about the invertibility of the operator in the 
left-hand side of Eq. ([?T[) is unnecessary. 



V. DIAGRAMMATIC REPRESENTATION 

An alternative method to obtain the spin-flip scatter- 
ing rate involves a direct calculation of the diagrams 
shown in Fig. [Tj These diagrams are superconducting 
cooperon propagators: they are built of the retarded 
Green function for particles [with dispersion ^(p)] and 
the retarded Green function for holes [with dispersion 
— ^(p)]. The latter can be converted into the advanced 
function for particles at the opposite energy, which we 
shall use below. These cooperon soft modes naturally 
arise in the diagrammatic expansion of the non-linear 
sigma model discussed in Sec. IIIIl Thus the diagrams of 
Fig. [T) directly correspond to the sigma-model diagrams 

of Fig. m 

The cooperon modes are massless in the absence of su- 
perconducting correlations (i.e., become singular if both 
the total external momentum and energy are zero). Once 
the spin flip is taken into account, the cooperon acquires 
a mass which is directly related to the spin-flip rate. 

Superconducting correlations in the system (e.g., due 
to the proximity effect) also produce a mass for the 
cooperon. We shall neglect this effect in the calculation 
of the spin-flip rate due to the conditions ([2|). All the 
mechanisms resulting in a mass of the cooperon, includ- 
ing the spin flip, are weak and their contributions may 
be calculated independently. 

In order to calculate the effective spin-flip rate, we shall 
evaluate the two diagrams depicted in Fig. [T) at the zero 
total momentum and zero energy. The diagrams (a) and 
(b) yield the following contributions to the cooperon self 
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energy: 



and 



(27r)3 



(27r)3 



G«(p + q)G«(p)G-4(-p) 



(59) 



respectively, where G^''*(p) — [— C(p) i ^]~^ ^^'^ 
the zero-energy retarded and advanced Green functions. 
C{q) denotes the cooperon containing only the potential 
impurities, 



Ciq) 



1 + B{q) + B^{q) 



2tivt 1 - B{q) 



where 



B{q) 



2ttvt J (27r)3 



G«(p + q)G-4(-p) 



(60) 



(61) 



is a single ladder rung containing a disorder line and two 
Green functions. The factor of 2 in Eq. ([55)1 comes from 
two possible diagrams of type (b) . 

The total spin- flip scattering rate is the sum of Eqs. 
([58|) and ([59|) with a coefficient that can be easily deter- 
mined from comparing to the limit of pointlike magnetic 
impurities.— In that limit, Fij{q) is actually independent 
of the momentum q and only the term (|58p contributes, 
yielding 



7^^^) - 47rVT2 I d\F,,{r). 



(62) 



By comparing with Eq. (jlip . we arrive at the result for 
the spin-flip rate: 



(63) 



Now the calculation of the spin-flip rate can be con- 
veniently performed in the two overlapping regimes: the 
medium-to- long-range magnetic correlations (a S> kp^) 
and the short-to-medium-range magnetic correlations 
(a < I). 

In the medium-to-long-range regime, the integrals in 
Eqs. 1(551) dSH) are restricted to q <^ kp and can be 
performed by using the integration over ^ in the vicinity 
of the Fermi surface. First, we use the identity 



G«(p)G^(-p) = [G«(p) - G^(-p)] 



(64) 



and discard all the integrals containing only retarded or 
only advanced Green functions (since they have all the 



poles lying in the same half-plane of the variable ^) . This 
allows us to re-express Eqs. ([55)1 and ((55)) in terms of the 
function B{q) defined in Eq. ((6T|) : 



7(a) 



,3 / 



— F,^{q)Biq), (65) 



1(b) 4^^^ / (2^)3 l-B{q) 



As a result, Eq. ((63|) leads to 

d^q F,^{q)B{q) 
(27r)3 1-S((7) 



(66) 



(67) 



An explicit calculation of B{q) (by using an integration 
over ^) gives 



B{q) 



1 



2'KVT 



d3p/(27r) 



(-^-vq+^)(-e-27) 

1 \ arctan(gZ) 



1 -I- iTvq/^ 



9/ 



(68) 



(the averaging here is over the Fermi surface). Substi- 
tuting this expression into Eq. (j67p , we obtain the result 
()15p for the medium- and long-range-correlation regimes. 

In the short-to-medium-range regime, the integrals 
(155)) and ((55| are determined by the momenta q 3> Z"'^, 
which allows us to neglect the contribution ()59p in favor 
of ((58| . If the correlation length of the magnetic disorder 
becomes comparable to kp^^ the integrals of G^G^ and 
G^G^ cannot be neglected any more, and the integral 
((55)) should be calculated in a different way. The main 
contribution to the p integral comes from the intersection 
of the two mass shells of the "width" l^^ shifted by the 
vector q. Using the inequality q ^ we approximate 
G^(p)G'^(— p) by the delta- function and obtain 

/ 7^G''(p + q)G^(-p-q)G«(p)G^(-p) 



d^p 

(27r)3 



S[^{p)] S[aP + q)] 



ql 



9{2kF - q) . (69) 



Substituting this expression into Eq. ((55)) . we arrive at 
the short-to-medium-range crossover result ()13p . 

Of course, in the quasiclassical limit kpl ^ 1 consid- 
ered in this paper, one is allowed to combine the two 
overlapping regimes into a single formula 

pij ^ [ <P<^ F,j{q) arctan(g/) 

J (2^)3 gZ-arctan(gO ' ^ ^ 

q<2kF 

which reproduces both the short-to-medium- and 
medium-to-long-range crossover results. 
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VI. CONCLUSIONS 

To summarize, we have analyzed the effect of magnetic 
inhomogeneities in disordered superconducting systems. 
Our approach covers magnetic disorder of various corre- 
lation lengths, and thus extends the theory of Abrikosov 
and Gor'kov for magnetic impurities,— as well as some 
earlier studies of superconductivity in systems with in- 
homogeneous mag netismJiii^i^ The main conclusion of 
our work is that if the correlation length of the mag- 
netic disorder is much shorter than all the macroscopic 
scales in the problem [condition ([2])], then the effect of 
the magnetic disorder may be incorporated as an effec- 
tive local spin-flip rate (in the same form as for magnetic 
impurities^ ) . We have obtained exact expressions for the 
effective spin-flip rate, under the assumption of a Gaus- 
sian magnetic disorder. 

While the exact expressions for the effective spin-flip 
rate are probably of mainly academic interest, we believe 
that our results will be helpful for estimating the spin-flip 
effects induced by inhomogeneities in various experimen- 
tal setups. As an example of such an application, we con- 
sider the experiments on SFS 7r-junctions, where the spin 
flip plays an important role^'^ (it manifests itself in the 
difference between the length scales involved in the de- 
cay and oscillations of the critical current as a function of 
the junction thickness). If we apply our estimates to the 
experimental data reported in Ref. d (assuming Sh ^ h), 
then we arrive at the estimate of the disorder correlation 
scale a ~ 2 nm. Note that this correlation scale is of the 
order of the length scales associated with the uniform 
component h of the magnetic field and with the result- 
ing spin-flip rate, thus, this example is at the borderline 



of applicability of our theory. The estimated size of in- 
homogeneities is apparently too small for domains (in 
recent experiments on CuNi films similar to those used 
in the 7r-j unctions, domains of size about 100 nm have 
been reported^). However, our estimates are consistent 
with earlier indications of clusters of magnetic Ni atoms 
in such alloys*^'^^ (inhomogeneities inside the domains). 

Finally, we would like to comment on comparison be- 
tween the effective spin- flip rates due to two cases of inho- 
mogeneous magnetization; disordered and periodic ones. 
The disordered case is considered in the present paper, 
while specific realizations of periodic magnetic structures 

were studied before in Refs. Ellilllilii. The obtained 

results for the spin-flip rate are all of the same order 
of magnitude, differing only by numerical factors. This 
suggests that there is probably no qualitative difference 
between the effective spin-flip rates in disordered and pe- 
riodic magnetic structures, as long as the characteristic 
length scale of inhomogeneities is sufficiently small. 
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